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Abstract 

The accurate solution of some of the main problems in numerical linear algebra 
(linear system solving, eigenvalue computation, singular value computation and 
the least squares problem) for a totally positive Bernstein- Vandermonde matrix is 
considered. Bernstein- Vandermonde matrices are a generalization of Vandermonde 
matrices arising when considering for the space of the algebraic polynomials of 
degree less than or equal to n the Bernstein basis, a widely used basis in Computer 
Aided Geometric Design, instead of the monomial basis. 

Our approach is based on the computation of the bidiagonal factorization of a 
totally positive Bernstein- Vandermonde matrix (or its inverse) by means of Neville 
elimination. The explicit expressions obtained for the determinants involved in the 
process makes the algorithm both fast and accurate. 

Key words: Bernstein- Vandermonde matrices, total positivity, Neville elimination, 
high relative accuracy. 



1 Introduction 



The design of accurate and efficient algorithms for structured matrices is a 
relevant field in numerical linear algebra which in recent years has received a 
growing attention (see, for instance, the recent survey paper [6] and references 
therein). An important condition for an algorithm to be accurate is the so 
called NIC (no inaccurate cancellation) condition, which is emphasized in [6]: 
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NIC: The algorithm only multiplies, divides, adds (resp., substracts) real 
numbers with like (resp., differing) signs, and otherwise only adds or substracts 
input data. 

In particular, in Section 2.3 of [6] several different classes of structured matrices 
are considered, among them the class of totally positive matrices. In our work 
we will consider one special class of totally positive matrices which we have 
recently analyzed in [T7] in the context of linear system solving: Bernstein- 
Vandermonde matrices. 

Bernstein- Vandermonde matrices are a generalization of Vandermonde ma- 
trices arising when considering for the space of the algebraic polynomials of 
degree less than or equal to n the Bernstein basis, a widely used basis in Com- 
puter Aided Geometric Design due to the good properties that it possesses 
(see, for instance, [^|7] ). The explicit conversion between the Bernstein basis 
and the power basis is exponentially ill-conditioned as the polynomial degree 
increases [8], and so it is very important that when designing algorithms for 
performing numerical computations with polynomials expressed in Bernstein 
form, all the intermediate operations are developed using this form only. 

As recalled in [6] , a crucial preliminary stage of the algorithms considered there 
for the class of totally positive matrices is the decomposition of the matrix 
as a product of bidiagonal factors, and that is our main contribution for the 
case of Bernstein- Vandermonde matrices: the fast and accurate computation 
of this bidiagonal factorization. 

Working along this line, in this work we present algorithms for linear sys- 
tem solving, eigenvalue computation, singular value computation and the 
least squares problem which are based on the bidiagonal factorization of the 
Bernstein- Vandermonde matrix or of its inverse. Factorizations in terms of 
bidiagonal matrices are very useful when working with Vandermonde [2][T3"] . 
Cauchy [3], Cauchy- Vandermonde [20f2~T] and generalized Vandermonde ma- 
trices [5]. 

As for the other stages of our algorithms, we will use several algorithms de- 
veloped by Plamen Koev [T1"]I15II16] . 

The rest of the paper is organized as follows. Some basic results on Neville 
elimination and total positivity are recalled in Section 2. In Section 3 the 
bidiagonal factorization of a Bernstein- Vandermonde matrix and of its inverse 
are considered. In Section 4 the algorithms for solving the four above problems 
are presented. Finally, Section 5 is devoted to illustrate the accuracy of the 
algorithms by means of some numerical experiments. 



2 



2 Basic facts on Neville elimination and total positivity 



To make this paper as self-contained as possible, we will briefly recall in this 
section some basic results on Neville elimination and total positivity which 
will be essential for obtaining the results presented in Section 3. Our notation 
follows the notation used in [9J and [10]. Given k, n E N (1 < k < n), Qk,n 
will denote the set of all increasing sequences of k positive integers less than 
or equal to n. 

Let A be an I x n real matrix. For k < I, m < n, and for any a G Qk,i and 
(3 G Qm,ni we will denote by A[a|/3] the sub matrix k x m of A containing the 
rows numbered by a and the columns numbered by (3. 

The fundamental tool for obtaining the results presented in this paper is the 
Neville elimination PfTU] . a procedure that makes zeros in a matrix adding to 
a given row an appropriate multiple of the previous one. We will describe the 
Neville elimination for a matrix A = {(kj)i<i<i;i<j<n where I > n. 

Let A = {di,j)i<i<i;i<j<n be a matrix where I > n. The Neville elimination of 
A consists of n — 1 steps resulting in a sequence of matrices A := A\ — > A2 — > 
. . . — > A n , where A t = {o^})i<i<i;i<j< n has zeros below its main diagonal in 
the t—1 first columns. The matrix A t+ i is obtained from A t (t = 1, . . . , ri) by 
using the following formula: 



a (m) - 

a i,3 



afj , if i < t 



<j - Kt/<4-iM-ij > if i>t + lmdj>t+l (2.1) 
, otherwise. 



In this process the element 

Pi,j ■= of) 1 < 3 < n, j <i <l 

is called pivot of the Neville elimination of A. The process would break 
down if any of the pivots pij (1 < j < n, j < i < I) is zero. In that case 
we can move the corresponding rows to the bottom and proceed with the new 
matrix, as described in [9]. The Neville elimination can be done without row 
exchanges if all the pivots are nonzero, as it will happen in our situation. The 
pivots are called diagonal pivots. If all the pivots pij are nonzero, then 
Pi,i = CLi,i Vi and, by Lemma 2.6 of [9J 

det A\i — j + 1, . . . , ill, . . . , ?1 , . 

Pi, = A , AV ~zi — ~rpr — tt K j < n, j < i < I. 2.2 

det A[i - j + 1, 1|1, ... ,j - 1} 
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The element 



'1,3 ~ 



Pi,3 



1 < 3 < n, j < i < I 



(2.3) 



is called multiplier of the Neville elimination of A. The matrix U := A n is 
upper triangular and has the diagonal pivots in its main diagonal. 

The complete Neville elimination of a matrix A consists on performing the 
Neville elimination of A for obtaining U and then continue with the Neville 
elimination of U T . The pivot (respectively, multiplier) of the complete 

Neville elimination of A is the pivot (respectively, multiplier) (j, i) of the 
Neville elimination of U T , if j >i. When no row exchanges are needed in the 
Neville elimination of A and U T , we say that the complete Neville elimination 
of A can be done without row and column exchanges, and in this case the 
multipliers of the complete Neville elimination of A are the multipliers of the 
Neville elimination of A if i > j and the multipliers of the Neville elimination 



A matrix is called totally positive (respectively, strictly totally positive) if all 
its minors are nonnegative (respectively, positive). The Neville elimination 
characterizes the strictly totally positive matrices as follows [9] : 

Theorem 2.1. A matrix is strictly totally positive if and only if its complete 
Neville elimination can be performed without row and column exchanges, 
the multipliers of the Neville elimination of A and A T are positive, and the 
diagonal pivots of the Neville elimination of A are positive. 

It is well known [4] that the Bernstein- Vandermonde matrix is a strictly totally 
positive matrix when the nodes satisfy < x\ < %2 < ■ ■ ■ < xi + i < 1, but this 
result is also a consequence of our Theorem 3.2. [19J. 



3 Bidiagonal factorizations 



The Bernstein basis of the space H n (x) of polynomials of degree less than or 
equal to n on the interval [0, 1] is: 



of A T if j > i. 




i = 








4 



The matrix 



.4 



(»>? 



is the (/ + 1) x (n + 1) Bernstein- Vandermonde matrix for the Bernstein basis 
B n and the nodes {%i}i<i<i+i- 



From now on, we will assume < X\ < x 2 < • • • < Xi + i < 1. In this situation 
the Bernstein- Vandermonde matrix is a strictly totally positive matrix [1] and 
the following two theorems hold: 



Theorem 3.1. Let A = ( 

^i,j)i<i,j<n+i be a Bernstein- Vandermonde matrix 
whose nodes satisfy < x\ < x 2 < ■ ■ ■ < x n < x n+ i < 1. Then A" 1 admits a 
factorization in the form 



A — G\G 2 ---G n D F n F n _x • • • Fi, 



(3-1) 



where Fi (i — 1, . . . , n) are (n + 1) x (n + 1) bidiagonal matrices of the form 



1 

1 



1 



— m 



i+2,i 



(3.2) 
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Gf (1 < i < n) are (n + 1) x (n + 1) bidiagonal matrices of the form 



G 



1 

1 



1 

-m i+2 ,i 1 



(3.3) 



y — m„ +lj j 1 , 

(i — 1, . . . , n), and D is a diagonal matrix of order n + 1 

D = diag{pi i i,p 2 ,2, • • • ,p n +i,n+i}- (3.4) 

rriij are the multipliers of the Neville elimination of the Bernstein- Vandermonde 
matrix A, and have the expression 



rriij 



(1 _ a;.)n-J+l(l - x-^.) ni=l(^ - Xi- k ) 



(1 - Xi_i)"-J'+ 2 ni= 2 (^-i - Xi_ fc ) 
where j = 1, . . . , n and i — j + 1, . . . , n + 1. 



(3.5) 



mjj are the multipliers of the Neville elimination of A T and their expression 
is 

(n — i + 2) • xj 



where j = 1, . . . , n and i = j + 1, . . . , n + 1. 



(3.6) 



Finally, the zth diagonal element of D is the diagonal pivot of the Neville 
elimination of A and its expression is 



ir fc -\(i - x fc ) 



(3.7) 



for z = 1, . . . , n + 1. 

Proof. It can be found in [17]. □ 

Theorem 3.2. Let A = (aij)i<i<.i+i-i<j< n +i be a Bernstein- Vandermonde 
matrix for the Bernstein basis E n whose nodes satisfy < x\ < x 2 < ■ ■ ■ < 
Xi < Xi + i < 1. Then A admits a factorization in the form 



A — FiFi-x ■ ■ ■ F\DG\ ■ ■ ■ G n -\G n 



(3.8) 
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where Fi (1 < % < I) are (I + 1) x (/ + 1) bidiagonal matrices of the form 



Fi 



1 

1 



V 



mti-i 1 



(3.9) 



Gf (1 < i < n) are (n + 1) x (n + 1) bidiagonal matrices of the form 



G 



1 

1 



1 



and D is the (I + 1) x (n + 1) diagonal matrix 

^ = (^i,j)l<i<«+l;l<i<n+l = diag{pi j i,P2,2) • • • ,Pn+l,n+l}- 



(3.10) 



(3.11) 



The expressions of the multipliers m it j (J = 1, . . . , n + 1; 2 = j + l,...,/ + l) 
of the Neville elimination of A, the multipliers m^j (j = 1, . . . , n; % = j + 
1, . . . ,n + 1) of the Neville elimination of A T , and the diagonal pivots P^i 
{i = 1, . . . , n+ 1) of the Neville elimination of A are also in this case the given 
by Eq. (3.5), Eq. (3.6) and Eq. (3.7), respectively. 

Proof. It can be found in IT9l. □ 



It must be observed that in the square case, the matrices Fj (i = 1, . . . , /) 
and the matrices Gj (j = 1, . . . ,n) that appear in the bidiagonal factoriza- 
tion of A are not the same bidiagonal matrices that appear in the bidiagonal 
factorization of A" 1 , nor their inverses (see Theorem 3.1 and Theorem 3.2). 
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The multipliers of the Neville elimination of A and A T give us the bidiagonal 
factorization of A and A' 1 , but obtaining the bidiagonal factorization of A 
from the bidiagonal factorization of A^ 1 (or vice versa) is not straightforward. 
The structure of the bidiagonal matrices that appear in both factorizations 
is not preserved by the inversion, that is, in general, F' 1 (i = 1, . . . ,/) and 
Gj 1 (j = l,...,n) are not bidiagonal matrices. See [TTJ for a more detailed 
explanation. 

A fast and accurate algorithm for computing the bidiagonal factorization of 
the totally positive Bernstein- Vandermonde matrix A and of its inverse (when 
it exists) has been developed by using the expressions (3.5), (3.6) and (3.7) 
for the computation of the multipliers m^j and TOjj, and the diagonal pivots 
Pi t i of its Neville elimination [r7]fT9] . Given the nodes {xj}i<j</ + i G (0, 1) and 
the degree n of the Bernstein basis, it returns a matrix M G R( z+1 ) x (™ +1 ) such 
that 





= Pi,i 


i = l,.. 


. ,n + 1, 






= rriij 


J = l,-- 


.,n + l; 


i = j + l,...,l + l, 




= rhj,i 


z = l,.. 


■,n; j = 


i + 1, . . . ,n + 1. 



The algorithm, which we have called TNBDBV, does not construct the Bernstein- 
Vandermonde matrix, it only works with the nodes {£i}i<i<z+i- Its computa- 
tional cost is of 0(ln) arithmetic operations, and has high relative accuracy 
because it only involves arithmetic operations that avoid inaccurate cancella- 
tion (see [T7] for the details). The implementation in Matlab of the algorithm 
in the square case can be taken from [14]. 

Remark. The algorithm TNBDBV computes the matrix M, denoted as BT>(A) 
in [T5], which represents the bidiagonal decomposition of A. But it is a remark- 
able fact that the same matrix BT>(A) also serves to represent the bidiagonal 
decomposition of A~ x . 



4 Accurate computations with Bernstein- Vandermonde matrices 

In this section four fundamental problems in numerical linear algebra (lin- 
ear system solving, eigenvalue computation, singular value computation and 
the least squares problem) are considered for the case of a totally positive 
Bernstein- Vandermonde matrix. The bidiagonal factorization of the Bernstein- 
Vandermonde matrix (or its inverse) let us to develop accurate and efficient 
algorithms for solving each one of these problems. 

Let us observe here that, of course, one could try to solve these problems by 
using standard algorithms. However the solution provided by them will gener- 
ally be less accurate since Bernstein- Vandermonde matrices are ill conditioned 
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(see [1 Tj ) and these algorithms can suffer from inaccurate cancellation, since 
they do not take into account the structure of the matrix, which is crucial in 
our approach. 

4-1 Linear system solving 

Let Ax = b be a linear system whose coefficient matrix A is a square Bernstein- 
Vandermonde matrix of order n + 1 generated by the nodes {£j}i<j< n +i, where 
< xi < . . . < x n+ i < 1. An application which involves the solution of this 
type of linear systems has been presented in [18J. 

The following algorithm solves Ax = b accurately with a computational cost 
of 0{n 2 ) arithmetic operations (see [I7J for the details): 

INPUT: The nodes {xi}i<i< n +i and the data vector b G R n+1 . 

OUTPUT: The solution vector x G R n+1 . 

- Step 1: Computation of the bidiagonal decomposition of A' 1 by using 
TNBDBV. 

- Step 2: Computation of 

x = A- l b = G 1 G 2 ■ ■ ■ G n D~ l F n F n _i ■■■F 1 b 

Step 2 can be carried out by using the algorithm TNSolve of P. Koev |14j . 
Given the bidiagonal factorization of the matrix A, TNSolve solves Ax = b 
accurately by using backward substitution. 

Several examples illustrating the good behaviour of our algorithm can be found 
in [UJ. 

4-2 Eigenvalue computation 

Let A be a square Bernstein- Vandermonde matrix of order n + 1 generated 
by the nodes {xj}i<j< n+ i, where < X\ < . . . < x n+ i < 1. The following 
algorithm computes accurately the eigenvalues of A. 

INPUT: The nodes {£j}i<j< n+ i. 

OUTPUT: A vector x G R n+1 containing the eigenvalues of A. 

- Step 1: Computation of the bidiagonal decomposition of A" 1 by using 
TNBDBV. 
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- Step 2: Given the result of Step 1, computation of the eigenvalues of A by 
using the algorithm TNEigenvalues. 

TNEigenvalues is an algorithm of P. Koev [15J which computes accurate eigen- 
values of a totally positive matrix starting from its bidiagonal factorization. 
The computational cost of TNEigenvalues is of 0(n 3 ) arithmetic operations 
(see |15j) and its implementation in Matlab can be taken from [13]. In this 
way, as the computational cost of Step 1 is of 0(n 2 ) arithmetic operations, 
the cost of the whole algorithm is of 0(n 3 ) arithmetic operations. 



4-3 The least squares problem 

Let A G R( /+1 ) x ( n+1 ) be a Bernstein- Vandermonde matrix generated by the 
nodes {xi}i<i<z+i, where < x\ < . . . < X[ + i < 1 and I > n. Let b G R/ +1 
be a data vector. The least squares problem associated to A and b consists on 
computing a vector x G R n+1 minimizing || Ax — b 1 1 s . 

Taking into account that in the situation we are considering A is a strictly 
totally positive matrix, it has full rank, and the method based on the QR 
decomposition due to Golub [12] is adequate [T]. For the sake of completeness, 
we include the following result (see Section 1.3.1 in [1]) which will be essential 
in the construction of our algorithm. 

Theorem 4.1. Let Ac = f a linear system where A G R( /+1 ) x(n+1 ), / > n, 
c G R n+1 and / G R z+1 . Assume that rank(A) — n + 1, and let the QR 
decomposition of A be given by 



where Q G R( /+1 ) x (' +1 ) is an orthogonal matrix and R G R( n+1 ) x ( n+1 ) i s a n 
upper triangular matrix with positive diagonal entries. Then the solution of 
the least squares problem min c || / — Ac H2 is obtained from 



where di G R n+1 , d 2 G R' n and r = / — Ac. In particular || r H2HI ^2 

The following algorithm, which is based on the previous theorem, solves in 
accurate and efficient way our least squares problem: 
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INPUT: The nodes {£j}i<i<«+i, the data vector / and the degree n of the 
Bernstein basis. 

OUTPUT: The vector x G R n+1 minimizing || Ax — b W2 and the minimun 
residual r = b — Ax. 

- Step 1: Computation of the bidiagonal factorization of A by means of 
TNBDBV. 

- Step 2: Given the result of Step 1, computation of the QR decomposition 
of A by using TNQR. 

- Step 3: Computation of 



The algorithm TNQR has been developed by P. Koev, and given the bidiagonal 
factorization of A, it computes the matrix Q and the bidiagonal factorization 
of the matrix R. Let us point out here that if A is strictly totally positive, then 
R is strictly totally positive. TNQR is based on Givens rotations, has a compu- 
tational cost of 0(l 2 n) arithmetic operations if the matrix Q is required, and 
its high relative accuracy comes from the avoidance of inaccurate cancellation 
|16j . Its implementation in Matlab can be obtained from |14j . 

As for the computational cost of the whole algorithm, it is led by the cost of 
computing the QR decomposition of A, and therefore, it is of 0(l 2 n) arithmetic 
operations. 

Some numerical experiments which show the good behaviour of our algorithm 
when solving problems of polynomial regression in the Bernstein basis have 
been presented in [T9j . 

4-4 Singular value computation 

Let A E R,( / + 1 ) x ( n+1 ) be a Bernstein- Vandermonde matrix generated by the 
nodes {£j}i<j<z+i, where < x± < . . . < xi + i < 1 and I > n. The following 
algorithm computes in an accurate and efficient way the singular values of A. 

INPUT: The nodes {xi}i<i<;+i and the degree n of the Bernstein basis. 




- Step 4- 

- Step 5: 



Solution of the upper triangular system Rc = d\. 
Computation of 
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OUTPUT: A vector x £ R n+1 containing the singular values of A. 

- Step 1: Computation of the bidiagonal decomposition of A by using TNBDBV. 

- Step 2: Given the result of Step 1, computation of the singular values by 
using TNSingularValues. 

TNSingularValues is an algorithm of P. Koev that computes accurate singular 
values of a totally positive matrix starting from its bidiagonal factorization 
|15j . Its computational cost is of 0(ln 2 ) and its implementation in Matlab 
can be found in [14] . Taking this complexity into account, the computational 
cost of our algorithm for computing the singular values of a totally positive 
Bernstein- Vandermonde matrix is of 0(ln 2 ) arithmetic operations. 



5 Numerical experiments 



In this last section we include two numerical experiments which illustrate the 
high relative accuracy of the algorithms we have presented for the problems 
of eigenvalue computation and singular value computation. Numerical experi- 
ments for the cases of linear system solving and of least squares problems have 
been included in IT71[T9l. 



Example 5.1. Let B20 be the Bernstein basis of the space of polynomials 
with degree less than or equal to 20 on [0, 1] and letA be the square Bernstein- 
Vandermonde matrix of order 21 generated by the nodes: 

22 ^ 20 18 16 14 ^ 12 ^ 10 8 6 4 ^ 2 ^ 42 ^ 38 34 ^ 30 

}^.<C — < — < — < — <- 
26 22 18 14 10 6' 

The condition number of A is ^{A) = 1.9e + 12. In Table 1 we present the 
eigenvalues Aj of A and the relative errors obtained when computing them by 
means of: 

(1) The algorithm presented in Section 4.2 (column labeled by MM). 

(2) The command eig from MATLAB. 

The relative error of each computed eigenvalue is obtained by using the eigen- 
values calculated in Maple 10 with 50-digit arithmetic. 

Example 5.2. Let i3 15 be the Bernstein basis of the space of polynomials 
with degree less than or equal to 16 on [0, 1] and let A £ 21x16 Bernstein- 
Vandermonde matrix generated by the same nodes as in Example 5.1. Its 
condition number is ^{A) = 5.3e + 08. In Table 2 we present the singular 
values Oi of A and the relative errors obtained when computing them by 
means of 
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MM 


eig 


l.Uc -f- uu 


n 
u 


A do 1 f; 

4.Ue — 10 


o.4e — Ul 


/I do 1 fi 

4.Ue — ID 


l.oc — 10 


9 8e ni 

Z.oe — Ul 


Q Qc 1 fi 

y.ye — id 


3 ne i ^ 

O.Uc — 10 


Z. 16 — Ul 


n 
u 


3 7c 1 ^ 
o. I c — 10 


l.Ze — Ul 


a ric 1 fi 
D.Uc — ID 


3 Re 1 ^ 

O.Oc — 10 


fi fie 09 
D.Dc — UZ 


8 ^e 1 fi 
O.Oc — ID 


X. ^c 1 ^ 

O.Oc — 10 


3 8c 09 
O.oe — UZ 


o. ( c — ID 


1 9o 1 /I 
l.Zc — 14 


9 9c D9 
Z.Zc — UZ 


A 8c 1 fi 
4.oe — 1U 


9 Dc 14 
Z.Uc — 14 


q /ic n3 
y.4e — uo 


n 
u 


/I fie 1/1 

4.0c — 14 


/i fie n3 
4.0e — uo 


1 Qe 1 fi 

i.ye — id 


Q 3e 1/1 

y.oe — 14 


1 ^e D3 
1.0c — Uo 


9 Qc 1 fi 

z.ye — id 


file 

0. 1c — 14 


0.y6 — U4 


1 8e 1 fi 

l.oc — ID 


9 ^e 13 

Z.Oc — lo 


1. 1 e — U4 


3 9o 1 fi 
O.Zc — ID 


1 3e 19 

l.oc — 1Z 


4.oe — UO 


Q Ac 1 fi 

y.4e — id 


3 fie 19 
O.Oc — 1Z 


i ^c 

l.oc — uo 


11c 1 ^ 
1.16 — 10 


fi Qe 19 

o.ye — iz 


i 7c nfi 
1. 1 e — uo 


file 1 fi 
0.1c — 10 


C 1„ 11 

0.1c — 11 


c a c n7 
0.0c — U f 


n 
u 


1 Qe in 

i.ye — iu 


3 ^e ns 

O.Oc — UO 


9 8c 1 ^ 

Z.Oc — lO 


1 8e DQ 
i.oe — uy 


Lie -08 


5.9e - 16 


5.9e - 09 


2.7e - 10 


2.1e - 15 


6.4e - 08 


1.3e - 12 


9.0e - 16 


1.0e-05 



Table 1 

Example 5.1: Eigenvalues of a Bernstein- Vandermonde matrix of order 21. 

(1) The algorithm presented in Section 4.4 (column labeled by MM). 

(2) The command svd from MATLAB. 

The relative error of each computed singular value is obtained by using the 
singular values calculated in Maple 10 with 50-digit arithmetic. 

Acknowledgements. This research has been partially supported by Spanish 
Research Grant MTM2006-03388 from the Spanish Ministerio de Education 
y Ciencia. 
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MM 


svd 


i p.o _i_ nn 

1.0c "T uu 


1 /1c 1 p, 

i.4te — id 


~\ Ac ~\P. 
1.4c — ID 


l.Zc "T UU 


Q Ac 1 fi 
y.4t6 — ID 


q Ac ^P. 
y.4tc — id 


^ So ni 
O.oc — Ul 


O.oc — ID 


O.oc — ID 


a So ni 

41. oc — Ul 


O.oc — ID 


9 ^o ~\P. 
Z.OC — ID 


9 ^o m 

Z.Oc — Ul 


9 9c 1 fi 
Z.Zc — ID 


n 
u 


i So ni 

l.oc — Ul 


r no 1 k 

D.Uc — ID 


A ^o 1 
4t.0c — ID 


r ^o no 
O.oc — uz 


9 9o 1 fi 
Z.Zc — ID 


n 
u 


^ ^o 09 
0.0c — UZ 


n 

U 


i no 1 ^ 

l.Uc — 10 


7 fie n 1 ? 

/ .0c — UO 


1 1 o 1 fi 

1 . 16 — ID 


Q So 1 ^ 

y.oc — 10 


q /!„ n^ 
0.41c — Uo 


1 1 o 1 fi 
1.1c — ID 


9 no 1 ^ 
z.uc — 10 


^ Q e nzL 
o.yc — U4t 


1 Ac 1 fi 
1.4tC — ID 


1 ^o 1 A 
1.0c — 14L 


i 'ic n/i 

l.OC — U4t 


9 no i r 

Z.Uc — ID 


8. Ac 1 ^ 
0.4tc — 10 


7 ^c nfi 

/ .Oc — Uu 


q rio i fi 

y.uc — iu 


1 fic 19 
l.Uc- — 1Z 


1.6e - 06 


1.32e - 16 


4.3e - 12 


4.8e - 08 


1.8e - 16 


3.9e - 10 


3.0e - 09 


2.9e - 15 


2.5e - 10 



Table 2 

Example 5.2: Singular values of a Bernstein- Vandermonde matrix 21 x 16. 
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